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We compute the contribution of classical fields to the second Fourier coefficient 
(772) of the azimuthal gluon distribution at large transverse momentum in heavy ion 
collisions. We find that the classical contribution to the flow alone cannot account 
for the experimentally observed behavior of V2{pt) at large transverse momentum pt. 



I. INTRODUCTION 



The azimuthal anisotropy of the particles produced in heavy ion collisions is quantified 
by the Elliptic flow parameter v2. It is defined to be the second Fourier moment of the 
azimuthal distribution, 



V, = (cos(20)) ^ — ;^ ■ (i) 



The RHIC (Relativistic Heavy Ion Collider) experiments have measured the centrality de- 
pendence of f 2 and find a value for V2 ^ that is significantly larger for non-central collisions 
than previous measurements at lower energies 0. Hydrodynamic models do a reasonable 
job of describing the centrality and mass dependence of V2 H, §]• The RHIC experiments 
have also measured V2{pt), which is defined to be 

r,#cos(20); '^'^ 

•f—TT dypt dpt dtj) 



^2 IP* J - TiT-Ti rf^TV ■ K"^) 



The data |l|, |] show that, for a wide range of centralities, V2{pt) rises up to 1.5 GeV, flattens 
and then remains flat up to the highest transverse momentum measured p^ ^ Q GeV. For 
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peripheral collisions (specifically 34%-85% central collisions), V2{pt) can be as large as 25% 
at Pi ~ 3 GeV. 

Hydrodynamic models describe the initial rise in V2{pt) but overpredict the data at large 
Pt. The models predict a continued rise except perhaps in a Blast wave parametrization 
with extreme assumptions 0. Partonic descriptions which account for the quenching of 
high Pt "jets" as they traverse the hot QCD matter do not describe the data either [0]-unless 
unusually large partonic cross-sections are assumed p. 

In this note, we discuss the computation of V2 in the Color Glass Condensate (CGC) 



approach to nuclear collisions 0, |I0|, [1^. At very high energies, the density of partons per 



unit area in the colliding nuclei becomes very large and saturates at a scale Qs ^ ^qcd [113 



For RHIC, estimates give Qs ~ 1 — 2 GeV [jT2|, |T^, |T6l. The typical occupation number of 
partons is l/as{Qs) > 1- Thus classical methods can be used to describe nuclear collisions at 
high energies. The CGC approach has been used to describe successfully various aspects of 



the bulk properties of charged hadrons in heavy ion collisions, including the multiplicity [12 



the centrality dependence of the multiplicity [0], the energy and rapidity dependence [|T7 
and the pt dependence of inclusive hadron spectra [0. 

Elliptic flow, as deflned in Eq. 1, is dominated by contributions from small momenta. 
In the classical approach, this involves computing diagrams to all orders in Qs/pt- This 
computation is very difficult to perform analytically. It can be computed numerically [ P0[ . 
However, when the transverse momenta is large, pt ^ Qs-, the classical equations can be 
linearized and the contribution to V2{pt) can then be computed analytically. This computa- 
tion is performed in this paper. We find that the contribution to V2{pt) from classical gluon 
flelds at large momenta is small and differs signiflcantly from the measured distribution at 
high Pt. The origin of the empirically observed behavior of V2{pt), at large pt, must arise 
from a source outside the classical perturbative expansion. Interestingly, even though the 
contribution from high momenta is small, this contribution is generated at very early times. 

The paper is organized as follows. In the following section, we set up the problem of 
computing the classical gluon distribution at large transverse momentum. In section 3, we 
perform the computation. A brief flnal section summarizes our results. Some details of the 
computation are relegated to an appendix. 
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II. THE CLASSICAL GLUON DISTRIBUTION 

Within the framework of the McLerran-Venugopalan model, Kovner, McLerran, and 
Weigert (KMW) (see also Refs. 0| and Jl^) computed the classical spectrum of hard 



gluons (with pt > Qs) radiated from a collision of two infinitely large nuclei with uniform 
color charge distributions. In order to calculate the azimuthal distribution of the radiated 
gluons, we generalize this work to a finite nucleus with a non-uniform color charge distribu- 
tion. Gradients in the charge distribution then give rise to space-momentum correlations of 
the radiated gluons. The typical momentum scale of the produced gluons is ~ Qs{xt). A 
typical spatial gradient is ~ 1/R. Therefore, the typical space-momentum correlation is of 
order ~ (QgR)^^. The final expression for t>2 is actually of order ~ (QgR)^^ since elliptic 
fiow responds to the quadrupole moment of the color charge distribution. 

We first review the perturbative solutions to the classical Yang Mills equations in order 
to establish the notation. In the McLerran-Venugopalan model, two nuclei (A and B) collide 
at very high energies. The valence partons generate classical fields at very small Bjorken 
X. The subsequent evolution of these classical fields is determined by the equations of 
motion. Therefore, the solution to the equations of motion is a functional of the valence 
color charge density (per unit area per unit rapidity) in nuclei A and B - pA^Xt) and psixt), 
respectively. The expectation of any quantity O is first expressed in terms of the classical 
solution and then averaged over all possible valence charge distributions with a weight 
function. In the McLerran-Venugopalan model, a Gaussian ansatz is taken for the weight 
function. Specifically, for a functional 0{pa,Pb) we have the following average: 

(0)^ = / [DpA][DpB]OipA,PB)exp(^- J d'x,^^^^^ - J . (3) 

At large transverse momentum, the density of charges is small and the field strength is 
weak. Therefore the classical Yang-Mills equations can be linearized and the solution can 
be expressed as a functional of the sources in closed form. The final distribution of gluons 



can then be related to the averages of the gluon field. The analysis of [|10| shows that the 
distribution of radiated gluons is given by 

dN 



dydPkt 



oc / £xtd'^yt(^''''''e-'^'y'i:i 



[Xt) 



d^A,d'(j)B (yt) 



t 



p 



x{6'^''^ + e'^e^^), (4) 
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where 0a (0b) is the 2-d Coulomb potential associated with the color charge distribution 

Pa (pb), 

- Vt4>a{xt) = gpA- (5) 

In what follows, we evaluate the average over p and and perform the necessary Fourier 
transforms in order to determine the gluon distribution to order ~ (Qs-R)"^- With this 
gluon distribution, we determine the azimuthal anisotropy, V2{pt)- 

III. COMPUTATION OF V2{k) FOR k > 

The averages over the color sources in Eq. |^ can be performed most simply in Fourier 
space. Below, we adopt the convention that repeated spatial indices are integrated over and 
repeated Fourier indices are integrated over with the appropriate two dimensional measure, 
/ (^^- In Fourier space, the average of 0{pa,Pb) is then 

{0)^ = J [Dp~A][DpB]0{pA,PB)exp{-TT{pA)l{CA)kkiPA)k'-Tr{pB)l{C^^^ (6) 

where {pA)k = e^^'^pA^x) is the Fourier transform of Pa{x) and {CA)kk' is given by 

(CaU' = e^-^e-^^'-' . (7) 
p\{x) 

Then pairwise correlations for nucleus A are given by 

((Pa)Upa).) = {Ct)u = e'^'-^>^pl{x) ^ p\ik - I) . (8) 
Fourier transforming each of the 0's, the gluon distribution of Eq. ^ becomes 

Contracting the pA variables and the pb variables we find 



oc 



A^^di + kYih + ky{CA%i,iCB')ih+km+k) . (10) 



dyd^k /f/i(/i + fc)2(/3 + fc)2 

Specializing now to the case where the target and projectile are identical, restoring the 
imphcit integration, and using {Cj^)ki = (C'^^)fc« = /^^(/c — /) we obtain 
dN f f d% d% ^^^ki^^ij^ki 



dyd^k J J {2TiY{2TiYllll{h + kY{h + k) 
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For the purposes of illustration, we consider a Gaussian model charge distribution 



Ax) 



111 RxRy 



(12) 
(13) 



where Rx and Ry denote the transverse spatial extent (in the x and y direction) of the 
overlap region of the two nuclei and is the net color charge squared in the overlap region. 



For this model distribution, we see that 



/i^(/i — I3) is a sharply falling function. Thus 



/i^(/i — /a) is approximately a delta function. Indeed, when ^"^{li — l 



is a delta function 



Eq. reduces to the result of ||TOl for the gluon distribution. We therefore approximate 
the integral over /i^ as 







1 (27r)2 





2 5.,5.J,=o 



(14) 



In writing the above equation we have assumed that the charged distribution 
symmetric in z and defined the moments of this distribution as 

(Pz 



f^^iz] 



IS 



(/^);f 

For the problem at hand, we have 

dN 



dycPk 



where f{k,z) is given by 



f{k,z)^ 



(27r) 



d'^z 
(2^ 



ij kl 



^2(^) 



(15) 
(16) 



^i\z) f{k,z) 



(17) 



-p*(p - z)''{p + ky{p - z + ky 



(27r)2 p2(p — zy(p + ky(p — z + ky 
We may write down the general form of the second derivatives of f{k, z) 



(18) 



' d^f{k,z) _ 6'^ d^f{k,zy 
^ dzidzj 2 d^^d^^ 

( dy{k,zy 



z=0 



2=0 



(19) 
(20) 
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Now the elliptic flow parameter f 2 is given by 



/gcos(20); '^^ 

2n dyk dk d<f> 



V2{y,k)^ ' ^;^^y^^''^^ . (21) 



Substituting Eq. |14| and Eq. |T9| into this deflnition of V2, we obtain to leading order in 

""'^y^^^- 8(/x4)(o) f{k,z = Q)- ^^^^ 
f{k, z = 0) and f^'^\k'^) are calculated in appendix A. The results are logarithmically diver- 
gent. As discussed in this is a consequence of the weak fleld expansion; higher order 
non-linear corrections cut off the behavior in the infra-red at a scale cts/x. To leading log 
accuracy we flnd 

^f^-^^")- MiJ'"^^^^^) '''' 

With these expressions we obtain our flnal result for f 2 

It is remarkable (see the appendix) that there are no power or even logarithmic divergences 



from the ratios of integrals in our expression for V2- For the Gaussian distribution of Eq. ^ 
the moments of the charge distribution are easily calculated. The resulting elliptic flow is 

The contribution of momenta k > Qg to the integrated V2 (Eq. |1|) is 

For realistic values of ~ 1 GeV for RHIC and i?^ ~ 3 fm for a peripheral collision, 
this contribution is a fraction of a percent. The signiflcant contribution of classical flelds 
to V2 comes from momenta k < Qg. Interestingly, our results suggest that V2 for k > Qg 
is generated very early. The spatial and temporal components decouple already at proper 
time r = 0. The temporal contributions are Bessel functions Jo{kT) and Ji{kT) which can 
be linearized into plane waves for times r < 1/Qs when k» Qs- 
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IV. SUMMARY 



In this paper, we have computed analytically the contribution to the azimuthal anisotropy 
{v2{k)) from classical fields at large transverse momenta, k"^ > Q^. We find that the az- 
imuthal anisotropy has the simple form shown in Eq. ^ for a Gaussian charge distribution. 
Clearly, this behavior disagrees with the data and another mechanism must be found to ex- 
plain the large momentum behavior. Non-fiow azimuthal correlations may provide a natural 
explanation [^. An interesting feature of our result is that the temporal dependence of 
the perturbative classical fields decouples already at proper time r = |T0[- Thus as Eq. ^ 
suggests, the momentum anisotropy is generated at very early times in the collision. This 
analysis may thus explain why (in numerical simulations of the Yang-Mills equations pO| ) 
a significant anisotropy is seen at early times after the collision. 



V. ACKNOWLEDGEMENTS 



We would like to thank D. Kharzeev, Y. Kovchegov, A. Krasnitz, L. McLerran, Y. Nara 
and K. Tuchin for useful discussions. R.V and D.T.'s work was supported by DOE Contract 
No. DE-AC02-98CH10886. R.V would also like to acknowledge RIKEN-BNL for support. 



APPENDIX A: EVALUATION OF INTEGRALS 



To complete the evaluation of f 2 we expand Eq. as a function of z and extract the 
coefficients of the symmetric traceless tensor structures. We find after some algebra 

= + + (A2) 

where the integrals are given by 

^^'^'^ = / (270^2(^+^)6 (^*^^- - f) (A3) 

~ J {27r)^ p^{k + pr[ 2 P 2) ^^^^ 
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f{k,z = 0) is logarithmically divergent at p = and p = —k. It is identical to the expression 



obtained previously in pO| , |T^. As discussed there and shown explicitly in Ref. |T^ this 
infra-red divergence is an artifact of the weak field expansion. To leading log accuracy we 
find 



Consider now the sum of the integrals in Eqs. While each of these individually have 
power divergences, taken together, the result is again only logarithmically divergent. Note, 
as previously, the integrals are divergent at p = and p = —k and one obtains logarithmic 
contributions from both regions to the final result. For illustration, consider the contribution 
to the sum from the divergence around p = 0. is logarithmically divergent, is 

finite, and appears power divergent. However, is actually only logarithmically 

divergent upon integration over the angle between p and k. Contracting both sides of Eq. 



with 2(^||i — ^) and keeping only logarithmically divergent terms we obtain the contribution 
of the p = divergence to and respectively to be 



k'k^ 




k^ 


~ T 


k'y 




A;2 


~ Y 



(27r)2p2(fc + p)6 (277) A;4 ^\{asiif, 



o(T\ — \ - ( 2p2cos(20pfc) _ 1 1 ( k^ \ 

^ [ 2 j,^,^^^^ - J (2vr)2 p^ik + py ^ (2.)fc4^"HK/^)V ^ ^ 

where (j)pk is the angle between p and k. A detailed analysis shows that the contribution to 

the sum + + from the logarithmic divergence around p = —k is identical 

to the sum of the divergent pieces in Eqs. \K7\ . 

Substituting the results of the above analysis into the expression for f'^'^^k'^) we find 

With these expressions for f{k,z = 0) and f^'^\k'^)^ the final result for V2 given in Eq. ^ 
follows. 
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